Dynamics and flow-induced phase sepa- 
ration in polymeric fluids 

Peter D. Olmsted 

The past few years have seen many advances in our un- 
Q"derstanding of the dynamics of polymeric fluids. These include 
(^improvements on the successful reptation theory; an emerging 
0"molecular theory of semifiexible chain dynamics; and an un- 
derstanding of how to calculate and classify "phase diagrams" 
£3br flow-induced transitions. Experimentalists have begun map- 
I ^ing out the phase behavior of wormlike micelles, a "living" 
polymeric system, in flow: these systems undergo transitions 
^jnto shear-thinning or shear-thickening phases, whose variety is 
remarkably rich and poorly understood. Polymeric ideas must 
r~ fo e extended to include the delicate charge and composition 
effects which conspire to stabilize the micelles and are strongly 
influenced by flow. 

-t-i . 

C&ddress 

£3)epartment of Physics and Astronomy & IRC in Polymer 
I S cience and Technology, University of Leeds, Leeds LS2 9JT, 
~3JK; p . d . olmstedgleeds .ac.uk 

O ' 

o : 

Introduction 



f-~Jhe field of non-linear rheology is roughly fifty years old and 
CNfar from mature. Non-linear fluids often display elastic ef- 
cts, and have effective viscosities which depend on stress 
strain rate. Polymeric liquids typically shear- thin [Jjj], al- 
inhough branched polymers thicken dramatically in extensional 
Q^low; colloidal suspensions of platelike particles (clays) typi- 
ng ally shear-thicken; and solutions of surfactant (soap) molecules 
C\&an shear-thicken or shear-thin. 

'"J" 1 , Non-linear rheology and polymer dynamics are immense 
'"^fields, and in this short (subjective) review I will focus on a 
C3ew subfields. I will discuss recent advances in modelling the 
Qlynamics of flexible and semifiexible polymer melts, including 
linear and complex topologies; and then review progress in our 
. ^p'owledge of the surfactant wormlike micelle system, to which 
^^oncepts from polymer dynamics have been successfully ap- 
plied. Unlike conventional polymers, the micellar microstruc- 
^ure can change qualitatively in flow conditions; these transi- 
tions have many features in common with equilibrium phase 
transitions, and have excited great interest. A good collection 
of results from a wide range of complex fluids may found in 



2 Polymer Melts and Solutions 

The most accepted molecular model for the dynamics of flex- 
ible entangled polymers has been the Doi-Edwards (DE) the- 
ory, based on de Gennes' reptation concept, in which polymers 
are envisaged to occupy "tubes" that model entanglement con- 
straints. This has done a reasonable job in predicting linear 



rheology, with a few notable exceptions such as the failure to 
predict the scaling of the zero frequency viscosity as M 3A , al- 
though recent work suggests that contour length [||| fluctuations 
are a key to this puzzle. 

Recent work in the linear regime includes applying the tube 
picture to molecules with complex branched topologies [H]. 
Star polymers afford stringent tests of tube model ideas, because 
diffusion is dominated by arm retraction in its tube, which is 
exponential in the retraction potential. Inclusion of higher or- 
der Rouse modes along with "dynamic dilution" of the tube has 
led to remarkably good agreement with experiment |Q . In ad- 
dition to star polymers, molecular models have recently been 
developed for progressively more complex topologies, paving 
the way for understanding the flow behavior of industrially- 
important long-chain-branched polymers, which strain-harden 
in extensional flow while softening in shear flow [|«4]|. 

Although the DE tube picture works well in the linear 
regime, it has several defects at high strain rates, particularly 
in steady shear. For example: experiments show a slightly in- 
creasing plateau shear and an increasing normal stress for strain 
rates 7 above an inverse reptation time r" 1 , while theory pre- 
dicts a decreasing stress a ~ 7 _1 and a constant normal stress; 
and the DE theory predicts a high strain rate viscosity which de- 
creases with molecular weight, while experiments merge onto a 
molecular weight-independent curve. The defect in DE theory 
is that, as the tube representing entanglement constraints rotates 
into the flow, the entrapped polymer feels a reduced stress and, 
at strain rates above the inverse tube relaxation time r" 1 , re- 
mains oriented and presents a decreasing stress with increasing 
strain rate. Although corrections due to tube stretching [||] have 
accounted for some problems in startup flows, this only applies 
near 7 ~ r" 1 and still predicts a pronounced stress maximum. 
Another mechanism is needed to relax the chain and hence in- 
crease the stress by providing more misaligned material for the 
flow to "grip". The key is believed to lie in "convected con- 
straint release", whereby the entanglement (tube) mesh convects 
away at high strain rates, leaving a relaxed coil. Early applica- 
tions of this idea [Q, have cleared up several problems with 
the DE theory, although the theory still predicts a slightly de- 
creasing stress with strain rate, so it should be regarded at pro- 
visional. 

While moderately successful molecular theories of flexible 
polymer and rigid rod dynamics have existed since the 70's, the 
study of semifiexible polymers (in which d < L p < L, where 
d is the polymer diameter, L p the persistence length, and L the 
length) is quite young, mainly due to the severe mathematical 
difficulties in treating the bend degrees of freedom and length 
constraint. However, with increasing attention being paid to bio- 
logical polymers such as actin [^|, a deeper understanding of the 
dynamics of semifiexible polymer solutions is emergi ng. D irect 
imaging of tagged fluorescent polymers is possible [»1C], and 
several techniques have been developed for measuring elastic 
moduli, including direct (torsional oscillator flTT|]) and indirect 
(from various optical techniques [»12, |»13 |). Experiments in- 
dicate the vestiges of a plateau modulus, less pronounced than 
that of flexible polymer solutions. For polymers with L p < L e , 
the distance between entanglements or confinement constraints 
(or deflection length [|l4|), one expects the behavior of flexi- 



ble solutions. However, for L p > L e one expects qualitatively 
different effects due to the perturbation of ben ding modes by 
tube constraints. While Odijk and Semenov Jl4], |j"5|] have stud- 
ied the dynamics and statistics of individual filaments, only re- 
cently have molecular theories for the stress response of entan- 
gled solutions emerged. Unlike flexible polymers, semifiexible 
chains have a bending energy which maintains L p , and one can 
distinguish between longitudinal and transverse conformational 
changes. Two pictures have emerged for the origin of elastic 
stress in concentrated solutions: Isambert and Maggs [ |l6| | ar- 
gued that semifiexible chains can slide along their tubes longi- 
tudinally, and relaxation only occurs when transverse motions 
allow escape from the tubes. MacKintosh et al. [17| argued 
that, if longitudinal motion is suppressed, then the modulus is 
due to the applied tension and the relaxation of bending modes 
(which are present in the quiescent state due to thermal fluctu- 
ations). In the case of solutions the former mechanism is ex- 
pected to hold at times longer than that on which chain tension 
can relax [ |T8| l . These pictures have been made more quantita- 
tive by Morse [••19|, who has developed a molecular theory at 



in quantitative agreement with experiment [»27]. 



the level of the Doi-Edwards theory and included the bending 
curvature explicitly in the expression for the microscopic stress 
tensor. 



Flow instabilities in Wormlike Mi- 
celles 



DE theory predicts a bulk flow instability in polymer melts 
which has not been seen; however, a suggestive instability 
known as the "spurt effect" has been seen in extrusion, in which 
the throughput increases dramatically above a critical pressure 
gradient, often accompanied by a spatial pattern in the extru- 
date [20 2TJ] . Current opinion is that this is a surface instabil- 
ity, although the picture is not settled [p2|]. However, there is 
a polymeric system which displays a well-documented bulk in- 
stability and has been the subject of intense investigation in the 
past decade. 



Certain aqueous surfactant solutions (e.g. cetylpyridinium 
chloride/sodium salicylate [CPCl/NASal]; cetyltrimethylammo- 
nium bromide(CTAB)/NaSal) self-assemble into flexible cylin- 
drical micelles with an annealed length distribution that can 
encompass polymeric dimensions (microns). These solutions 
comprise a surfactant (e.g. CPC1) and an ionizing salt (e.g. 
NaSal) which together determine micellar dimensions, flexibil- 
ity, and interactions. Salt and concentration effects are quite 
delicate, with Coulomb interactions playing an important and 
poorly-understood role. Micelle reaction kinetics introduce ad- 
ditional timescales to the Rouse and reptation times of conven- 
tional polymers: in the limit of fast breaking times the stress 
relaxation of entangled micelles often obeys a simple single ex- 
ponential ("Maxwell fluid"), and properties can be calculated 
quite confidently [ [23[ |24j|, in good agreement with experiment 
J25|]. Some non-linear properties can be calculated in this limit, 
and a maximum in the shear stress (analogous the stress max- 
imum in DE theory) is predicted at an inverse relaxation time 
(the geometric mean of the reptation and breaking times) [E6L 
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Figure 1: Model constitutive relation cr(j) for a shear-thinning 
complex fluid with a bulk instability (after Ref. [|29"||). The thick 
curve shows the kind of curve envisioned for wormlike micellar 
solutions: this one is from the Johnson-Segalmann model [63]. 
Portions of the curve with negative slope are mechanically un- 
stable. Both the Doi-Edwards theory for polymer melts and the 
Cates theory for micelles predict a constitutive curve qualita- 
tively like the thin curve, with a continuously decreasing stress 
above a critical strain rate j c . Convected constraint release [ |»8| | 
has reduced the instability in the polymer melt theory to a slight 
bump, while micellar theories still rely on an implicit solvent 
viscosity to stabilize the high strain rate branch (if there is in- 
deed one). The portion of the curve with negative slope is un- 
stable. Experiments show that, for applied strain rates 7 >j p , 
micellar solutions can phase separate under shear, attaining a 
unique steady state stress 17^(7) in both phases, with a portion 
of the sample at a high shear rate 72 and a portion at the low 
shear rate 71 = j p . A dependence of a p on 7 (not in this case) 
implies different compositions in the coexisting phases, so that 
different constitutive curves are connected; and all properties 



change with mean strain rate [»«36 



The non-linear rheology of micellar systems was first stud- 



ied by Rehage and Hoffman [»27], who discovered dramatic 



shear thinning (analogous to the DE instability) and shear- 
thickening, depending on the salt/surfactant/water composition. 
The shear-thinning systems were the first to be systematically 
studied (see Figure 1). Above a critical strain rate j p an ap- 
parent phase separation into macroscopic coexisting regions oc- 
curs, at a reproducible and history-independent stress a p , in 
which the high strain rate material is well aligned (typically 
birefringent) and the low strain rate material remains relatively 
disordered. The underlying flow curve has the stress maximum 
Omax (predicted by Cates [ |26| ] for semi-dilute systems), while 
the composite steady state flow curve has a plateau beginning at 
c p < Umax- [It is important to note that micellar systems have 
slow dynamics, and one can trap metastable states for a > <j p 
p9]].] This occurs in semi-dilute systems, of order a few per- 
cent surfactant pjj|], or in more concentrated systems (of order 
30%) with a nearby equilibrium nematic transition [[30| ^TJ|. In 
the former case the dynamic instability is believed to be poly- 



meric in nature [|26j], while the latter may be due to nematic 
effects (probably both effects are present). No theories exist 
for nematic transitions under shear in micelles, although recent 
work include s phas e diagrams for model rigid-rod suspensions 
in shear flow [ »»32 ], for which only a few results exist [Ref. [^| 
reported a shear-induced nematic transition in a liquid crystal 
polymer melt] . 

Shear banding can be inferred from rheological measure- 
ments and directly observed optically. Quantitative measure- 
ments include the fraction of material and degree of alignment 
in the two phases, inferred from neutron scattering [»34]; and 
the velocity profile, measured directly using magnetic resonance 
imaging [»35|. Shear-banding can incorporate different concen- 
trations in the two phases, which is expected when flow modifies 
intermicellar interactions (as near a nematic transition) rather 
than simply the micellar conformation (as might be expected 
in more dilute systems). A signature of this is a slope in th e 
"plateau" stress with increasing mean strain rate [ »»36 , »»32 |, 
indeed seen in concentrated solutions which often have an un- 
derlying nematic transition [|37l ^8 , 39, »34|. 



Groups have begun investigating metastability. Berret et 
al. @ examined slow transients in 10-20% CPCl/NaSal 
solutions. After increasing the strain rate into the two-phase 
region the stress decayed slowly in time from the underly- 
ing constitutive curve onto the stress plateau, with behavior 
(7 ~ cxp — {t/r (7)}° , {a = 2), which they interpreted as one- 
dimensional nucleation and growth. Grand et al. [»42| stud- 
ied transients in more dilute (~ 1%) CPCl/NaSal solutions and 
found similar stress decays, with a ~ (2, 2.5, 3), and r(-f) di- 
verging above or below (depending on composition) the strain 
rate 7 P at the onset of banding. They also performed controlled 
stress experiments, and discovered a stress aj ump > a p , below 
which the system remained on the low strain rate branch indefi- 
nitely and above which the system eventually jumped to the high 
strain rate branch. Their data suggested that some compositions 
behave "spinodal-like" and others behave "nucleation-like" (as 
Berret's did), but it is to early to completely embrace the lan- 
guage of first-order transitions (given, e.g., <Jj ump , which has 
no equilibrium analog). 

Fischer and Rehage showed how shear-thinning systems can 
be tuned by changing surfactant and salt composition, from a 
shear banding material to a material with a stress plateau (the 
rheological signature of banding) but without banding [»43|. 
The shear and normal stresses apparently follow the Giesekus 
model, which is one of the simplest non-linear constitutive equa- 
tions (comprising a Maxwell model with the simplest stress- 
dependent relaxation time). A molecular understanding for this 
behavior is lacking. 

With critical micelle concentrations of order a few parts per 
million, micelles entangle at astonishingly low dilutions. Amaz- 
ingly, systems which shear-thin at concentrations of a few per- 
cent can undergo a shear-thickening transition at fractions of 
< 0.1% [|»27l 144], |»»46| ]. This this shear-induced structure (SIS) 
is still undetermined; early suggestions for the mechanism in- 
cluded runaway micellar growth due to flow-alignment 
but the observed strain rate (of order inverse milliseconds) is 
much slower than the necessary micellar reorientation time, of 
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Figure 2: Qualitative constitutive curves for a shear-thickening 
system, such as the micellar system reported by Boltenhagen 
et al. 
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The thin curves denote single phase rela- 
tions (micellar liquid or a gel-like material); the thick curve de- 
notes the composite flow curve along which the material phase 
separates for controlled stress a > a c , and complete conver- 
sion to gel occurs for a = a max . Upon controlling the strain 
rate the system traverses between the two branches at j c . Note 
that, if compositions of the coexisting phases were identical, the 
controlled-stress composite flow curve would be vertical. 



order /is. It is probable that charge, which controls the dramatic 
increase in micellar length for concentrations near the overlap 
concentration p9[], plays an important role. Like the shear- 
thinning systems, macroscopic "phase separation" occurs; Hu 
et al. 



••46 1 found a gel-like phase that forms upon increasing 
the applied stress, with the mean strain rate decreasing (see Fig- 
ure 2) as more material turns into gel, and increasing again after 
complete conversion. The gel is observed to fracture in flow, 
and slightly shear-thins. Applying a strain rate above the critical 
strain rate induces immediate complete conversion. Attempts to 
visualize the SIS using cryo-TEM have given few clues to the 
microstructure |^|. Note that coexistence in the shear-thinning 
micelles occur under controlled strain rate conditions, while co- 
existence in this thickening system occurs for controlled stress; 
in both cases banding occurs in the radial direction, indicating 
banding at a common shear stress. These differences may be co- 
incidences of the constitutive behaviors of the coexisting phases, 
or due to whether stress or strain rate ultimately determines the 
SIS. 



Berret et al. [ 



•51]] studied cethyltrimethylammonium tosy- 
late (CTAT) micelles, and found shear-thickening phase sepa- 
ration under controlled strain rate conditions above a critical 
strain rate j c ~ tp - 55 (an increase in j c with </> was also found 
by Hu et al. [ |»»46 ]); this concentation dependence remains 
unexplained. The composite curve 0-^,(7 ) has a positive slope, 
in contrast to the S curve of Ref. [ |»«46 |, possibly because the 
"thick" phase is not thick enough; alternatively, phase separa- 
tion along the vorticity direction (at a common strain rate) would 
also be consistent with a positive slope da p /d r y for the com- 
posite flow curve [••32]. The SIS is shear-thinning, displays 
an oriented structure in neutron scattering, and does not have 
the extremely long recovery times found in Ref [••46|. Qual- 



itatively similar data were reported for a CTAB-NaTOS solu- 
tion by Harmann and Cressely [^21 We finish our (incomplete) 
zoo of micellar thickening transitions by mentioning yet another 
study on CPCl/NaSal: revisiting early experiments by Rehage, 
Wheeler et al. [ 44 1 found spatio-temporal instabilities consist- 
ing of dark and light oscillating vertical bands (in cylindrical 
Couette flow); these accompany formation and destruction of 
new microstructure (evident from turbidity), and are suggestive 
of a Taylor-Couette elastic instability [p3]. 

Athough wormlike micelles can have much simpler rheol- 
ogy than their polymer cousins due to the frequent presence of 
a single relaxation time in the fast breaking limit, this simplicity 
is delicate, and strong flows can dramatically affect the micelle 
microstructure. We are far from a general molecular theory for 
these transitions, and do not even know the nature of the mi- 
crostate in most cases. Continuum constitutive models may pro- 
vide some insight, although when these models su ccee d we usu- 
ally do not know why (e.g. The Giesekus model [»43]). A pop- 
ular constitutive model is the local Johnson-Segalman model 
[p4]l, which is relatively simple and displays the non-monotonic 
flow curve characteristic of shear-thinning micelles. Numerical 
calculations [ |55[ |56| ] resemble some startup experiments [p4|, 
and authors have attempted to determine plateau stress a p for 
the onset of banding in this model [|4[ |56]]. It has become appar- 
ent that this, or any local, model does not give a unique selected 
banding stress [p7[], and an additional assumption is necessary 



[p8, »»32]. There is growing consensus that non-local (i.e. gra- 



dient terms) contributions to constitutive equations supply an 
unambiguous determination of the plateau stress [••32]. 



4 Outlook 



Despite progress in understanding of flexible and semiflexible 
polymer dynamics, there is no shortage of problems for the 
immediate future. We lack a credible complete molecular 
understanding for any of these micellar flow-induced transti- 
tions: such a model must presumably include charge and 
concentration, as well as polymeric effects, to account for the 
(still unknown, in most cases!) structural changes under flow. 
By comparing and contrasting living and non-living polymers 
we may be able to extract important physics. Besides these 
systems, many other complex fluid systems undergo a variety of 
flow-induced phase transitions, and it seems reasonable to hope 
that this variety of transitions may be put on a common ground, 
akin to the thermodynamics of equilibrium phase transitions. 

I am indebted to J-F Berret, ME Cates, SL Keller, CYD 
Lu, FC MacKintosh, TCB McLeish, DJ Pine, and G Porte for 
much discussion and advice. 
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